To APPEAR IN THE ASTROPHYSICAL JOURNAL 

Preprint typeset using I^'T^]X style eniulatcapj v. 11/10/09 



O 
(N 



(N 

< 

6 



> 

o 
\o 
o 
vn 

i> 

o 
o 



% 



THE RELATION BETWEEN GAS AND DUST IN THE TAURUS MOLECULAR CLOUD 

Jorge L. Pineda^, Paul F. Goldsmith^, Nicholas Chapman^, Ronald L. Snell^, Di Li^, Laurent Cambresy'', 

Chris Brunt* 

ijet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109-8099, USA 
^Department of Astronomy, LGRT 619, University of Massachusetts, 710 North Pleasant Street, Amherst, MA 01003, USA 
•^Observatoire Astronomique de Strasbourg, 67000 Strasbourg, France 
^Astrophysics Group, School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL, UK 

To appear in the Astrophysical Journal 

ABSTRACT 

We report a study of the relation between dust and gas over a 100 deg^ area in the Taurus molecular 
cloud. We compare the H2 column density derived from dust extinction with the CO column density 
derived from the ^^CO and ^"^CO J = 1 — ?► lines. We derive the visual extinction from reddening 
determined from 2MASS data. The comparison is done at an angular size of 200", corresponding to 
0.14 pc at a distance of 140 pc. We find that the relation between visual extinction Ay and A^(CO) is 
linear between Ay ~ 3 and 10 mag in the region associated with the B213-L1495 filament. In other 
regions the linear relation is fiattened for Ay > 4 mag. We find that the presence of temperature 
gradients in the molecular gas affects the determination of A^(CO) by ~30-70% with the largest 
difference occurring at large column densities. Adding a correction for this effect and accounting 
for the observed relation between the column density of CO and CO2 ices and Ay, we find a linear 
relationship between the column of carbon monoxide and dust for observed visual extinctions up to 
the maximum value in our data ~ 23 mag. We have used these data to study a sample of dense cores 
in Taurus. Fitting an analytical column density profile to these cores we derive an average volume 
density of about 1.4 x 10** cm~^ and a CO depletion age of about 4.2 x 10^ years. At visual extinctions 
smaller than ~3mag, we find that the CO fractional abundance is reduced by up to two orders of 
magnitude. The data show a large scatter suggesting a range of physical conditions of the gas. We 
estimate the H2 mass of Taurus to be about 1.5 x 10^ M©, independently derived from the Ay and 
A^(CO) maps. We derive a CO integrated intensity to H2 conversion factor of about 2. 1x10^° cm~^(K 
km s^^)~^, which applies even in the region where the [C0]/[H2] ratio is reduced by up to two orders of 
magnitude. The distribution of column densities in our Taurus maps resembles a log-normal function 
but shows tails at large and low column densities. The length scale at which the high-column density 
tail starts to be noticeable is about 0.4 pc. 
Subject headings: ISM: molecules — ISM: structure 



1. INTRODUCTION 

Interstellar dust and gas provide the primary tools for 
tracing the structure and determining the mass of ex- 
tended clouds as well as more compact, dense regions 
within which new stars form. The most fundamental 
measure of the amount material in molecular clouds is the 
number of H2 molecules along the line of sight averaged 
over an area defined by the resolution of the observa- 
tions, the H2 column density, A^(H2). Unfortunately, H2 
has no transitions that can be excited under the typical 
conditions of molecular clouds, and therefore it cannot 
be directly observed in such regions. We have to rely on 
indirect methods to determine A^(H2). Two of the most 
common methods are observations of CO emission and 
dust extinction. 

Carbon monoxide (CO) is the second most abundant 
molecular species (after H2) in the Universe. Observa- 
tions of ^^CO and ^^CO together with the assumption of 
local thermodynamic equilibrium (LTE) and moderate 
^■^CO optical depths allow us to determine iV(CO) and, 
assuming an [C0]/[H2] abundance ratio, we can obtain 
A^(H2). This method is, however, limited by the sen- 
sitivity of the ^^CO observations and therefore is only 
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able t o trace large column densities. iGoldsmith et al.l 
(|200l used a 100 square degree map of ^^CO and "CO 
in the Taurus molecular cloud to derive the distribution 
of iV(CO) and iV(H2). By binning the CO data by exci- 
tation temperature, they were able to estimate the CO 
column densities in individual pixels where ^^CO but not 
^■^CO was detected. The pixels where neither ^^CO or 
^■^CO were detected were binned together to estimate the 
average column density in this portion of the cloud. 

Extensive work has been done to assess the reliability 
of CO as a tracer o f the column of H2 molecules (e.g. 
iFrerking et all Il982t iLanger et aLlUQSOt ). It has been 
found that iV(CO) is not linearly correlated with A^(H2), 
as the former quantity is sensitive to chemi cal effects such 
as CO depletion at high volume densities (iKramer et ahl 
[1991 ICaselh et al.l fl99l iTafalla et al.l 12002^ and the 
competition between CO f ormation and d estructi on at 
low-column densi ties (e.g. Ivan Dishoeck fc Black 119881 : 
iVisser et al.ll2009( ). Moreover, temper ature gradients are 
likely present in molecular clouds (e.g. lEvans et al.ll2001| ) 
affecting the correction of A^(CO) for optical depth ef- 
fects. 

The H2 column density can be independently in- 
ferred by measuring the optical or near-infrared light 
from background stars that has been extincted by the 
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dust present in the molecular cloud (jLada et alj 119941 : 
iCambresvl [1991 iDobashi et alll2005[ ). This method is 
often regarded as one of the most reliable because it 
does not depend strongly on the physical conditions of 
the dust. But this method is not without some uncer- 
tainty. Variations in the total to selective extinction and 
dust-to-gas ratio, particularly in denser clouds like those 
in Taurus, may introduce some uncertainty in the con- 
version of the infrare d extinction to gas column density 
(jWhittet et al.l[2001[ ). Dust emission has been also used 
to derive the column density of H2 (Langcr ct al. 1989). 
It is, however, strongly dependent on the dust temper- 
ature along the line of sight, which is not always well 
characterized and difficult to determine. Neither method 
provides information about the kinematics of the gas. 

It is therefore of interest to compare column density 
maps derived from ^-^CO and ^'^CO observations with 
dust extinction maps. This will allow us to character- 
ize the impact of chemistry and saturation effects in the 
derivation of A^(CO) and A^(H2) while testing theoreti- 
cal predictions of the physical processes that cause these 
effects. 

As mentioned before, CO is frozen onto dust grains in 
regions of relat ively low temperatu r e and larger volume 
densities (e.g. iKramer et all Il999t iTafalla et all 120021: 
iBergin et al.ll2002l). In dense cores, the co lumn densities 
of C ^ ^O j Bcrgin ct al.' '200?) a nd C^^O ('Kramer_elaLl 
119991 lAlve s ct al. 1999; Kainula inen et al . 2006) are ob- 
served to be linearly correlated with Ay up to ^10 mag. 
For larger visual extinctions this relation is flattened with 
the column density of these species being lower than that 
expected for a constant abundance relative to H2 . These 
authors showed that the C^^O and C^^O emission is op- 
tically thin even at visual extinctions larger than 10 mag 
and therefore the flattening of the relation between their 
column density and Ay is not due to optical depths ef- 
fects but to depletion of CO onto dust grains. These 
observations suggest drops in the relative abundance of 
C^^O averaged along the line-of-sight of up to a factor of 
~3 for visual extinctions between 10 and 30 mag. A sim- 
ilar result has been obtained from direct determinations 
of the column density of CO-ices based on abso rption 
studi es toward embedded and field stars (Chia r et al.l 
|1995[) . At the center of dense cores, the [C0]/][H2] ratio 
is exp ected to be reduced by up to five orders of magni- 
tude (IBergin fc Langeil ll997\ This has been confirmed 
by the comparison between observations and radiative 
transfer calculations of dust continu um and C^^O emis- 
sion in a sample o f cores in Taurus (jCaselli et al.l 119991 : 
ITafalla et al.ll20d^ . The amount of depletion is not only 
dependent on the temperature and density of the gas, 
but is also dependent on the timescale. Thus, determin- 
ing the amount of depletion in a large sample of cores 
distributed in a large area is important because it al- 
low us to determine the chemical age of the entire Tau- 
rus molecular cloud while establishing the existence of 
any systematic spatial variation that can be a result of a 
large-scale dynamical process that lead to its formation. 



At low column densities {A\ 
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3 mag. 



A^(CO) <W 



the relative abundance of CO 



and its isotopes are affected by the relative rates 
of formation and destruction, carbon isotope ex- 
change and isotope selective photodissociation by 



far-ultraviolet (FUV) photons. These effects can reduce 
[C0]/[H2] by up to t hree o rders of magnitude (e.g. 



van Dishoeck fc Blackl Il988t iLiszti l2007t IVisser et al.] 
2009() . This column density regime has been studied 



in dozens of lines-of-sight using UV a nd optical ab- 
sorption (e.g. JFederman et al.: 1980: Sheffer et al.l [20021: 
iSonnentrucker et al.l l2003t iBurgh et al.l |2007[) as well 
as in absor ption toward mm-wave continuum sources 
(iLiszt fc Lu cas 1991). The statistical method presented 
bv [Goldsmith et alj (|2008[ ) allows the determination 
of CO column densities in several hundred thousand 
positions in the periphery of the Taurus molecular cloud 
with A^(CO) ~ 10^'* — 10-^^cm~^. A comparison with 
the visual extinction will provide a coherent picture of 
the relation between A^(CO) and iV(H2) from diffuse to 
dense gas in Taurus. These results can be compared 
with theoretical predictions that provide constraints in 
physical parameters such as the strength of the FUV 
radiation field, etc. 

Accounting for the various mechanisms affecting the 
[C0]/[Il2] relative abundance allows the determination 
of the H2 column density that can be compared with 
that derived from Ay in the Taurus molecular cloud. 
It has also been suggested that the total molecular 
mass can be determined using only the integrated in- 
tensity of the ^^CO J = 1 ^- line together with the 
empirically-derived CO-to-Il2 conversion factor {Xqo = 
N(R2)/Igo)- The Xqo factor is thought to be depen- 
dent on the physical conditions of the CO-emitting gas 
([Malonev fc Black 1988) but it has been found to attain 
the canonical value for our Galaxy even in diffuse re- 
gions where the [C0]/][H2] ratio is st rongly affec ted by 
CO formation/destruction processes ([Lisztl [20071) . The 
large-scale maps of A^ (CO) and Ay also allow us to assess 
whether there is H2 gas that is not traced by CO. This 
so-called "dark gas" is suggested to account for a sub- 
stantial fraction of t he total molecular gas in our Galaxy 
([Grenier et al.l[2005[) . 

Numerical simulations have shown that the proba- 
bility density function (PDF) of volume densities in 
molecular clouds can be fitted b y a log-normal dis- 
tribution (e.g. lOstriker et al.ll2001l: INordlund fc Padoaiil 



\T9§^ ILi et all l2004HKlessenl 12000^ The shape of 
the distribution is expected to be log-normal as mul- 
tiplicative effects determine the volume dens i ty of 
a molecul ar cloud (Passot fc V azguez-Semadenil 119981 : 
Vazquez- Semadeni fc Garcia 200l|). A log-normal func- 
tion can also describe the distr ibution of column 
densities in a molecular cloud (jOstriker et al.l [20011 : 
[Vazquez-Semadeni fc Garcia[|200lD . For some molecular 
clouds the column density distribution can be well fitted 
by a log-normal (e.g. Wong et al. 2008; Goodman et al.l 
[20091 ). A study by Kainul ainen et al.[ ([2009D . however, 
showed that in a larger sample of molecular complexes 
the column density distribution shows tails at low and 
large column densities. The presence of tails at large col- 
umn densities seems to be linked to active star-formation 
in clouds. The Ay and CO maps can be used to deter- 
mine the distribution of column densities at large scales 
while allowing us to study variations in its shape in re- 
gions with different star-formation activity within Tau- 
rus. 

In this paper, we compare the CO column den- 
sity derived using the ^^CO and ^^CO data from 
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Right Ascension (J3000) 



Figure 1. Mask regions defined in the Taurus Molecular Cloud. 
Mask 2 is shown in black, Mask 1 in dark gray, and Mask in light 
gray. We also show th e 156 stellar members of Taurus compiled by 
ILuhman et al.l 120061) as white circles. 



iNaravanan et all 120081 (see also iGoldsmith et al.l l2008| l 
with a dust extinction map of the Taurus molecular 
cloud. The paper is organized as follows: In Section [5] 
we describe the derivation of the CO column density in 
pixels where both ^^CO and ^'^CO were detected, where 
^^CO but not ^'^CO was detected, as well as in the region 
where no line was detected in each individual pixel. In 
Section |3] we make pixel-by-pixel comparisons between 
the derived A^(CO) and the visual extinction for the large 
and low column density regimes. In Section [4. II we com- 
pare the total mass of Taurus derived from A^(CO) and 
Ay. We also study how good the ^^CO luminosity to- 
gether with a CO-to-H2 conversion factor can determine 
the total mass of a molecular cloud. We study the distri- 
bution of column densities in Taurus in Section 14.21 We 
present a summary of our results in Section [5] 

2. THE N(H2) MAP DERIVED FROM i^co AND i^co 

In the following we derive the column density of 
CO using the FCRA O 14-m ^^CO and ^=^C obser- 
vations presented by INaravanan et al.l 120081 (see also 
IGoldsmith et al.ll2008D . In this paper we use data cor- 
rected for error beam p ick-up using the method pre- 
sented by iBensch et al.l ()2001[ ). The correction proce- 
dure is described in Appendix [K\ The correction for er- 
ror beam pick-up improves the calibration by 25-30%. 
We also improved the determination of jWCO) compared 
to that presented by IGoldsmith et al.l (|2008| ) by includ- 
ing an updated value of the spontaneous decay rate and 
using an exact numerical rather than approximate an- 
alytical calculation of the partition function. The val- 
ues of the CO column density are about ^20 % larger 
tha n those presented bv IGoldsmith et aTl ()2008[ ) . Follow- 
ing IGoldsmith et al.l (|2008l ). we define Mask 2 as pixels 
where both ^^CO and ^^CO are detected, Maskl as pix- 
els where ^^CO is detected but ^^CO is not, and MaskO 
as pixels where neither ^^CO nor ^■^CO are detected. We 
consider a line to be detected in a pixel when its intensity, 
integrated over the velocity range between 0-12 km s~^, 
is at least 3.5 times larger than the rms noise over the 
same velocity interval. We show the mask regions in Fig- 
ure [TJ The map mean rms noise over this velocity range 



is CTT,;, =0.53Kkms-i for ^^CO and crT_.^^=0.23Kkms-i 

for ^■^CO. The map mean signal-to-noise ratio is 9 for 
i^CO and 7.5 for ^^CO. Note that these values differ 
slightly from those presented by Goldsmith et al. (2008) , 
as the correction for error beam pick-up produces small 
changes in the noise properties of the data. 

2.1. CO Column Density in Mask 2 
2.1.1. The Antenna Temperature 

When we observe a given direction in the sky, the an- 
tenna temperature we measure is proportional to the 
convolution of the brightness of the sky with the nor- 
malized power pattern of the antenna. Deconvolving 
the measured set of antenna temperatures is relatively 
difficult, computationally expensive, and in consequence 
rarely done. The simplest approximation that is made 
is that the observed antenna temperature is that coming 
from a source of some arbitrary size, generally that of the 
main beam, or else a larger region. It is assumed that 
the measured antenna temperature can be corrected for 
the complex antenna response pattern and its coupling 
to the (potentially nonuniform) source by an efficiency, 
characterizing the coupling to the source. This is often 
taken to be T^mb, the coupling to an uniform source of 
size which just fills the main lobe of the antenn a pat- 
tern. This was the approach used by Goldsmith et al.l 
(2008). In Appendix Elwe discuss an improved technique 
which corrects for the error pattern of the telescope in 
the Fourier space. This technique introduces a "corrected 
main-beam temperature scale", T,„b,c- We can write the 
main-beam corrected temperature as 



Tinb.c — Tq 



1 



;To/T„ _ I 



^To/Tbg 



1 



1-e- 



(1) 



where Tq — hv /k, Tex is the excitation temperature of the 
transition, Tf,g is the background radiation temperature, 
and T is the optical depth. This equation applies to a 
given frequency of the spectral line, or equivalently, to a 
given velocity, and the optical depth is that appropriate 
for the frequency or velocity observed. 

If we assume that the excitation temperature is inde- 
pendent of velocity (which is equivalent to an assumption 
about the uniformity of the excitation along the line-of- 
sight) and integrate over velocity we obtain 



rmb,c(w)dw 



ToC(Tex) 



=To/T„ 



- /"(l-e-^(^) 



)dv , (2) 



where we have included explicitly the dependence of the 
corrected main-beam temperature and the optical depth 
on velocity. The function C(Tox), which is equal to unity 
in the limit Ti,g — > 0, is given by 



c(re,) = 1 - 



„To/r„ 



1 



0T0/T5, 



(3) 



2.1.2. 



Ttie Optical Depth 

The optical depth is determined by the difference in the 
populations of the upper and lower levels of the transition 
observed. If we assume that the line-of-sight is charac- 
terized by upper and lower level column densities, Njj 
and Nl, respectively, the optical depth is given by 
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Figure 2. Correction factors for the relation between integrated 
main-beam temperature and upper level column density for a 
Gaussian velocity distribution of the optical depth. The dotted 
(blue) curve shows the correction factor obtained usi ng integrals of 
functions of the optical depth as given by Equation llSI The solid 
(red) curve shows the correction fa ctor employing the peak values 
of the functions, given by Equation llGI 



hvo 



T = -^^{ly) [NlBlu - NuBul] 
c 



(4) 



where vq is the frequency of the transition, (?!)(i^) is the 
hne profile function, and the B's are the Einstein B- 
coefficients. The line profile function is a function of the 
frequency and describes the relative number of molecules 
at each frequency (determined by relative Doppler veloc- 
ity). It is normalized such that J (t){i')di' — 1. For a 
Gaussian line profile, the line profile function at line cen- 
ter is given approximately by 0(fo) = I/'^^'fwhm, where 
SiypwHM is the full width at the half maximum of the line 
profile. 

We have assumed that the excitation temperature is 
uniform along the line of sight. Thus, we can define the 
excitation temperature in terms of the upper and lower 
level column densities, and we can write 



Ny _ gu_ 
Nl~ ql' 



-To /To, 



(5) 



where the g's are the statistical weights of the two levels. 
The relationship between the B coefficients, 



lets us write 



guBuL — QlBlu, 
hvoBuL(i>{v)Nu 



To/n, 



1 



(6) 



(7) 



Substituting the relationship between the A and B coef- 
ficients, 



A-UL = BjjL - 



Snhi 



(8) 




9 10 11 

^^CO Excitation Temperature 

Figure 3. Parameters of Mask 1 binned by ^^CO excitation tem- 
perature Tex. The bottom panel shows the derived H2 density {left- 
hand scale, squares) and the number of pixels in each Tex bin (right- 
hand scale, triangles). The most common Tex values are between 
5 — 9K. The middle panel shows the observed (left-hand scale, 
squares) and derived (right-hand scale, triangles) ^■^CO/^'^CO ra- 
tio. The top panel shows the derived ^-^CO column density as- 
suming a line width of lKms~^. The H2 density, ^^CO column 
density and derived ^^CO/^^CO ratio increase monotonically as a 
function of ^^CO excitation temperature. 



gives US 



'M 



c^AuL4>iv)Nu 



.To/T„ 



1 



(9) 



If we integrate both sides of this equation over a range 
of frequencies encompassing the entire spectral line of 
interest, we find 



T(i')dv 



c^AulN, 



UL^^U 



,To/T„ 



1 



(10) 



2.1.3. 



Upper Level Column Density 

It is generally more convenient to describe the optical 
depth in terms of the velocity offset relative to that of 
the nominal line center. The incremental frequency and 
velocity are related through dv = {c/uQJdu, and hence 
j T{v)dv = [c/vq) j T{v)dv. Thus we obtain 



T{v)dv 

We can rewrite this as 
1 



c'AulNu 



Sttz^o 



pTo/Te, 



1 



c^AulNu 



dTo/T(3x 



1 



STTt'o^ J T{v)dv 



(11) 



(12) 
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Table 1 

^^CO Excitation Temperature Bins in Mask 1 and Best Estimates of Their Characteristics 



J ex 

(K) 



i2cO/"CO 
Observed 



Number 
of Pixels 



n(H2) 
(cm^'^) 



N(CO)/Sv 
(10^6 cm-2/kms-i) 



Abundance Ratio 



5.5... 
6.5... 
7.5... 
8.5... 
9.5... 
10.5., 
11.5., 
12.5., 



21.21 


118567 


250 


0.56 


17.29 


218220 


275 


0.95 


14.04 


252399 


275 


1.6 


12.43 


223632 


300 


2.3 


11.76 


142525 


300 


3.6 


11.44 


68091 


400 


4.1 


11.20 


24608 


500 


5.3 


11.09 


6852 


700 


6.7 



30 
30 
30 
32 

40 
45 
55 



Substituting this into Equation ([2]), we can write an ex- 
pression for the upper level column density as 



Nn = 



STrki/^ 



hc^AuLC{T,^) 



J T(v)dv 



/(I 



^i^'))dv 



(13) 

For the calculation of the ^■^CO column densities (Sec- 
tion I2.1.4P we use a value for the Einstein A-coefiicient 
oiAuL =6.33x10"^ s-i (IGoorvitchlll99l . 

In the limit of optically thin emission for which t[v) 
■^ 1 for all V, and neglecting the background term in 
Equation pjPI. the expression in square brackets is unity 
and we regain the much simpler expression 



Nu{thin) 



SnkvQ 



h(?A 



UL 



T'mb,c(w)'iw 



(14) 



We will, however, use the general form of Nu given in 
Equation (fT3|) for the determination of the CO column 
density. 

We note that the factor in square brackets in Equa- 
tion (jl3|) involves the integrals of functions of the optical 
depth over velocity, not just the functions themselves. 
There is a difference, which is shown in Figure [21 where 
we plot the two functions 



CF{integral) 



J T{v)dv 



J {1 ^ e-^(^))dv ' 



and 



CF(peak) = - — 



To 



-TO ' 



(15) 



(16) 



as a function of the peak optical depth tq. There is a 
substantial difference at high optical depth, which re- 
flects the fact that the line center has the highest optical 
depth so that using this value rather than the integral 
tends to overestimate the correction factor. 



2.1.4. Total 



' CO column densities derived from ^^ CO and 
CO observations. 



In LTE, the column density of the upper level (J 
is related to the total ^^CO column density by 



N: 



"CO 



= iVi 



hBo.J{,I + l) 



U 



(2J-K1) 



(17) 



^ This usually does not result in a significant error since in LTE 
even in dark clouds Tox is close to 10 K as compared to Tf,g = 2.7 
K. Since Ti,g is significantly less than To, the background term is 
far from the Rayleigh-Jeans limit further reducing its magnitude 
relative to that of the first term. 



where _Bo is the rotational constant of ^■^CO (Bo = 5.51 x 
10^*^ s^^) and Z is the partition function which is given 

by 



,7=0 



(2J-t-l)e- 



-hBo(.7 + l) 



(18) 



The pa rtition function can be evaluated explicitly as a 
sum, but iPen ziasI (|1975D pointed out that for tempera- 
tures T 3> hBo/K , the partition function can be approx- 
imated by a definite integral, which has value kT/hBo. 
This form for the partition function of a rigid rotor 
molecule is almost universally employed, but it does con- 
tribute a small error at the relatively low temperatures of 
dark clouds. Specifically, the integral approximation al- 
ways yields a value of Z which is smaller than the correct 
value. Calculating Z explicitly shows that this quantity 
is underestimated by a factor of ~1.1 in the range be- 
tween 8 K to lOK. Note that to evaluate Equation (|18|) 
we assume LTE (i.e. constant excitation temperature) 
which might not hold for high-J transitions. The error 
due to this approximation is, however, very small. For 
example, for rox=10K, only 7% of the populated states 
is at J = 3 or higher. 

We can calculate the column density of ^'^CO from 
Equation (J17p determining the excitation temperature 
Tox and the "CO optical depth from i^co and "CO 
observations. To estimate Tex we assume that the ^^CO 
line is optically thick (r ^ 1) in Equation ([1]). This 
results in 



T,- 



5.53 



In 1 



5.53 



+0.83 



(19) 



where T^ ^ is the peak corrected main-beam bright- 
ness temperature of ^^CO. The excitation temperature 
in Mask 2 ranges from 4 to 19 K with a mean value of 
9.7 K and standard deviation of 1.2 K. 

Also from Equation ([T]), the optical depth as a function 
of velocity of the ^'^CO J = 1 — ?► line is obtained from 
the main-beam brightness temperature using 



r"(u) = -ln 



1- 



5.29 



j5.29/T„, _ ^ 



-0.16 



(20) 
where T^ ^ is the peak corrected main-beam brightness 

temperature of "'^^CO. We use this expression in Equa- 
tion (fT5)) to determine opacity correction factor. We 
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evaluate the integrals in Equation (jl5|) numerically. The 
correction factor ranges from 1 to ~4 with a mean value 
of 1.3 and standard deviation of 0.2. The ^^CO column 
density is transformed to ^^CO col umn density a ssuming 
a i^cO/i^CO isotope ratio of 69 CWilsonI 119991) . which 
should apply for the well-shielded material in Mask 2. 

2.1.5. Correction for Temperature Gradients along the Line 

of Sight 

In the derivation of the CO column density and its 
opacity correction we made the assumption that the gas 
is isothermal. But observations suggest the existence of 
core- to-edge tem perature differences in molecular clouds 
(e.g. lEvans et~a l. 2001) which can be found even in re- 
gions of only moderate radiation field intensity. There- 
fore the presence of temperature gradients might affect 
our opacity correction. 

We used the radiative tran sfer code RATRAN 
(jHogerheiide fc van der Ta3l2000[ ) to study the effects of 
temperature gradients on the determination of iV(CO). 
The modeling is described in the Appendix [Cl We found 
that using ^^CO to determine the excitation tempera- 
ture of the CO gas gives the correct temperature only 
at low column densities while the temperature is over- 
estimated for larger column densities. This produces an 
underestimate of the ^"^CO opacity which in turn affects 
the opacity correction of A^(CO). This results in an un- 
derestimation of A^( CO) . We derived a correction for 
this effect (Equation |C2| ') which is applied to the data. 

2.2. CO Column Density in Mask 1 

The column density of CO in molecular clouds is com- 
monly determined from observations of ^^CO and ^^CO 
with the assumption of Local Thermodynamic Equilib- 
rium (LTE), as discussed in the previous section. The 
lower limit of iV(CO) that can be determined is there- 
fore set by the detection limit of the ^^CO J = 1 — > 
line. For large maps, however, it is possible to determine 
A^(CO) in regions where only -'^^CO is detected in indi- 
vi dual pixels by using the statistical approach presented 
by [Goldsmith et al.l ()2008[ ) . In the following we use this 
approach to determine the column density of CO in Mask 
1. 

We compute the excitation temperature from the ^^CO 
peak intensities for all positions in Mask 1 assuming that 
the emission is optically thick. The Mask 1 data is then 
binned by excitation temperature (in IK bins), and the 
^■^CO data for all positions within each bin averaged to- 
gether. In all bins we get a very significant detection 
of "'^'^CO from the bin average. Thus, we have the exci- 
tation temperature and the observed ratio of integrated 
intensities (^^CO/^^CO) in each 1 K bin. Since positions 
in Mask 1 are distributed in the periphery of high ex- 
tinction regions, it is reasonable to assume that the gas 
volume density in this region is modest, and thus LTE 
does not necessarily apply, as thermalization would imply 
an unreasonably low gas temperature at the cloud edges. 
We therefore assume that ^^CO is sub-thermally excited 
and that the gas has a kin etic temperature of 15 K . We 
use the RADEX program (jvan der Tak et al.l [20071) . us- 
ing the LVG approximation, and the collision cross sec- 
tions from th e Leiden Atomic and Molecular Database 
(LAMDA: .Schoier et al.l [20051 ) . to compute line intensi- 



ties. The free parameters in the modeling are tempera- 
ture (T), density (n), CO column density per unit line 
width {N{C0)/5v), and the i^CO/^CO abundance ratio 
(i?). Since the excitation is determined by both density 
and the amount of trapping (N/Sv), there is a family of 
n — N{C0)/6v parameters that give the same excitation 
temperature. The other information we have is the ^^CO 
integrated intensity for the average spectrum in each bin. 
Thus the choice of n, N{CO)/Sv and R must reproduce 
the excitation temperature and the observed ^^CO/^'^CO 
ratio. Solutions also must have an optical depth in the 
^^CO J = 1 -^ Oofat least 3, to be consistent with 
the assumption that this isotopolog ue is optically thick. 
This is the same method used in Gol dsmith et al.l p008) . 
although this time we used the RADEX program and the 
updated cross-sections from LAMDA. 

In fact, at low excitation temperature the data can only 
be fit if the CO is strongly fractionated. At high excita- 
tion temperature we believe that the CO is unlikely to be 
fractionated, and thus, R must vary with excitation tem- 
perature. We chose solutions for Mask 1 that produced 
both a monotonically decreasing R with decreasing ex- 
citation temperature and a smoothly decreasing column 
density with decreasing excitation temperature. The so- 
lutions are given in Table 1 and shown Figure [31 The 
uncertainty resulting from the assumption of a fixed ki- 
netic temperature and from choo sing the best valu e for R 
is about a factor of 2 in A^(CO) (|Goldsmith et al.|[2008) . 

To obtain A^(CO) per unit line width for a given value 
of the excitation temperature we have used a non-linear 
fit to the data, and obtained the fitted function: 



/ N{CO) \ f_Sv_ 
\ cni^^ / \kms^^ 



6.5 X 10 



13 



-^ ex 

1^ 



(21) 



Wc multiply by the observed FWHM line width to de- 
termine the total CO column density. The upper panel 
in Figure [3] shows N{CO)/Sv as a function of Tex- 

2.3. CO Column Density in Mask 

To determine the carbon monoxide column density in 
regions where neither ^^CO nor ^-^00 were detected, we 
average nearly 10^ spectra to obtain a single ^^CO and 
^■^CO spectra. From the averaged spectra we obtain a 
^■^CO/^"^CO integrated intensity ratio of ~17. We need 
a relatively low R to reproduce such a low observed value. 
Values of _R = 25 or larger cannot reproduce the observed 
isotopic ratio and still produce ^^CO emission below the 
detection threshold. Choosing _R = 20 and a gas kinetic 
temperature of 15 K, we fit the observed ratio with n — 
100cm^3 andiV(CO) = 3 x W^ cm-"^ . This gives rise to 
a ^■^CO intensity of 0.7 K, below the detection threshold, 
however much stronger than the Mask average of only 
0.18 K. Thus, much of Mask must not contribute to 
the CO emission. In fact, only 26% of the Mask area 
can have the properties summarized above, producing 
significant CO emission. Therefore, the average column 
densit5l3 throughout Mask is 7.8 x 10" cm'^. 

^ Note th at the estima te of the CO column density in MaskO 
by Goldsmit h et al.l ll200g) did not include the ~26% filling factor 
we derived here and in consequence overestimated the CO column 
density in this region. 
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Figure 4. Maps of the CO column density {upper panel) and visual extinction {lower panel) in the Taurus Molecular cloud. The gray-scale 
in the A'^(CO) and Av maps is expressed as the square root of the CO column density and of the visual extinction, respectively. The angular 
resolution of the data in the figure is 40" for N{CO) and 200" for Ay- 
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log(N(CO)/cnT'^) 

Figure 5. Histogram of ttie ^^CO column density distributions 
in the Mask 0, 1, and 2 regions mapped in Taurus. The Mask 
is indicated by a vertical line at Af(CO) = 3 X 10^^ cm~'^ which 
represents the column density i n the CO-emitting region (26% of 
the area of Mask 0; see Section 12.31 1. Note that we have not yet 
corrected N{CO) in Mask 2 for the effect of temperature gradients 
in the opacity correction. 

Another option is to model the average spectra of ^^CO 
and ^^CO matching both the ratio and intensity. Since 
now, our goal is to produce CO emission with inten- 
sity 0.18 K, both i^CO and "CO wih be optically thin. 
Therefore we need an R that is equal to the observed 
ratio. For i? = 18, a solution with n = 100 cm~^, 
Sv = Ikms-i, and iV(CO) = T.SxlOi'^cm-^ fits both 
the ^^CO and ^'^CO average spectra for Mask 0. Note 
that this is very similar to the average solution (with a 
slightly larger R) that assumes that ~26% of the area 
has column density 3x10-'^^ cm~^ and the rest 0. Thus 
for a density of 100 cm^'^, the average CO column den- 
sity must be about 7.8x10-'^'* cm~^ in either model. Of 
course, if we picked a different density we would get a 
slightly different column density. As mentioned above, 
the uncertainty is iV(CO) is about a factor of 2. 

Note that the effective area of CO emission is uniformly 
spread over MaskO. We subdivided the ^^CO data cube 
in the MaskO region in an uniform grid with each bin 
containing about 10"* pixels. After averaging the spectra 
in each bin we find significant "'^^CO emission in 95% of 
them. 

3. COMPARISON BETWEEN Av AND N{^'^CO) 

In order to test our estimate of A^(CO) and assess 
whether it is a good tracer of A^(H2), we compare Mask 1 
and 2 in our CO column density map of Taurus with a 
dust extinction map derived from 2MASS stellar colors. 
Maps of these quantities are shown in Figure H) We also 
show in Figure [S] a histogram of the ^^CO column den- 
sity distributions in the Mask 0, 1, and 2 regions mapped 
in Taurus. The derivation of the dust extinction map is 
described in Appendix [B] The resolution of the map is 
200" (0.14pc at a distance of 140 pc) with a pixel spacing 
of 100". For the comparison, we have convolved and re- 
gridded the CO column density map in order to match 
this resolution and pixel spacing. 

3.1. Large N(^'^CO) Column Densities 

We show in Figure [6] a pixel-by-pixel comparison be- 
tween visual extinction and ^^CO column density. The 




10000 



E 



Figure 6. Comparison between the visual extinction derived 
from 2MASS stellar colors and the ^^CO column density derived 
from ^"^CO and ^^CO observations in Taurus. The dark blue line 
represents the ^^CO column density d erived from Ay as suming 
Af(H2)/Av = 9.4 X lO^'^cm-^mag-i IIBohlin et al.l \T97B) and a 
[CO]/[H2] abundance ratio of 1.1 X 10~^. The gray scale repre- 
sents the number of pixels of a given value in the parameter space 
and is logarithmic in the number of pixels. The red contours are 
2,10,100, and 1000 pixels. Each pixel has a size of 100" or O.OTpc 
at a distance of 140 pc. 

visual extinction and A'^(CO) are linearly correlated up 
to about Ay ~ 10 mag. For larger visual extinctions 
TV (CO) is largely uncorrelated with the value of Ay. In 
the range 3 < Ay < 10 mag, for a given value of ^v, 
the mean value of A^(CO) is roughly that expected for 
a [C0]/[H2] relative abundance of ^10~^ which is ex- 
pected for shielded regions (Solomon &: Klempereri[T972l : 
IHerbst fc KlemDereHll973D . Some pixels, however, have 
CO column densities that suggest a relative abundance 
that is reduced by up to a factor of ~3. In the plot 
we show lines defining regions containing pixels with 
Ay > 10 mag and with 3 < Ay < 10 mag and N{CO) > 
9 X 10^^ cm^^. In Figure[7]we show the spatial distribu- 
tion of these pixels in iV(CO) maps of the B213-L1457, 
Heiles's cloud 2, and B18-L1536 regions. White con- 
tours correspond to the pixels with Ay > 10 mag and 
black contours to pixels with 3 < Ay < 10 mag and 
iV(CO) > 9 X lO^^cm-^. Regions with Ay > 10 mag 
are compact and they likely correspond to the center of 
dense cores. The largest values of A''(CO), however, are 
not always spatially correlated with such regions. We no- 
tice that large iV(CO) in the Ay = 3 — 10 mag range are 
mostly located in the B213-L1457 filament. We study 
the relation between Ay and A^(CO) in this filament by 
applying a mask to isolate this region (see marked re- 
gion in Figure [7]). We show the relation between Ay 
and iV(CO) in the B213-L1457 filament in the left hand 
panel of Figure [5] We also show this relation for the 
entire Taurus molecular cloud excluding this filament in 
the right hand panel. Visual extinction and CO column 
density are linearly correlated in the B213-L1457 fila- 
ment with the exception of a few pixels that are located 
in dense cores (Cores 3, 6 and 7 in Table 2). Without 
the filament the N{CO)/Ay relation is linear only up to 
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Figure 7. N{CO) maps of the B213-L1495 (top), Heiles's cloud 
2 (middle), and B18-L1536 {bottom) regions. The white contours 
denote regions with Ay > 10 mag, while the black contours denote 
regions with Av < 10 mag and A'"(CO) > 9 X lO^'^cm"^ (see 
Figure |6]l. The blue contour outlines approximately the B213— 
L1457 filament. 



~4 magnitudes of extinction. In Section 13.1.11 we will 
see that the deviation from a linear N{CO)/Ay relation 
is mostly due to depletion of CO molecules onto dust 
grains. Depletion starts to be noticeable for Ay > 4 mag. 
Therefore, pixels on the B213-L1457 filament appear to 
show no signatures of depletion. This can be due either 
to the filament being chemically young in contrast with 
the rest of Taurus, or to the volume densities being low 
enough that desorption processes dominate over those of 
adsorption. If the latter case applies, and assuming a 
volume density of n(Il2) — 10"^ cm^"^ (low enough to not 
show significant CO depletion but still larger than the 
critical density of the "'^'^CO J = 1 — >■ line), this fila- 
ment would need to be extended along the line-of-sight 
by 0.9-3 pc for 3 < Ay < 10 mag. This length is much 
larger than the projected thickness of the B213-L1495 fil- 
ament of ~0.2pc but comparable to its length of ~7pc. 
We will study the nature of this filament in a separate 
paper. 

Considering only regions with Ay < 10 mag and 
A^(CO) > 10" cm-2 (see Section [S^) we fit a straight 
line to the data in Figure |6] to derive the [C0]/[Il2] 
relative abundance in Mask 2. A least squares fit rc- 
suhs in iV(C0)/cm-2 = (l.Ql ± 0.008) x lO^Mv/mag. 
Assuming that all hydrogen is in molecular form 
we can write the ratio betwe en H2 column den sity 
and color excess observed by iBohlin et al.l (|1978[ ) as 
A^(Il2)/EB-v=2.9xl02icm-2 mag~^ We combine this 
relation with the ratio of total to selective ex tinction 
i?v = Ay/EB-y c:i 3.1 (e.g. iWhitteti [20031) to ob- 
tain N{R2)/Ay = 9.4 X 102022cm-2 mag-^ Combin- 
ing the N (112) /Ay relation with our fit to the data, 
we obtain a [C0]/[Il2] relative abundance of 1.1x10^^. 
Note that, as discussed in Appendix [BJ grain growth 
would increase the value o f Ry up to ^4.5 in dense re- 
gions (jWhittet et al.|[200ll) . Due to this effect, we esti- 
mate that the derived Ay would increase up to 20% for 
Ay <10mag. This would reduce the N{}l2)/Ay con- 
version but also increase the Ay/N{CO) ratio. Thus the 
derived [C0]/[Il2] abundance is not significantly affected. 

3.1.1. CO depletion 

The flattening of the Ay-N{CO) relation for Ay > 
10 mag could be due to CO depletion onto dust grains. 
This is supported by obser vations of the pre-stellar core 
B68 by iBergin et al.l ()2002l ) which show a linear increase 
in the optically thin C-'^^O and C^^O intensity as a func- 
tion of Ay up to ~7mag, after which the there is a 
turnover in the intensity of these molecules. This is simi- 
lar to what we see in Figure[Bl Note, however. Ay alone is 
not the sole parameter determining CO freeze-out, since 
this process also depen ds on density and timescale (e.g. 
IBergin fc Langeilll997n. 

Following 'W hittet et al.l (|2010( ). we test the possibility 
that effects of CO depletion are present in our observa- 
tions of the Taurus molecular cloud by accounting for 
the column of CO observed to be in the form of ice on 
the dust grains. iWhittet et al.l ([^07) measured the col- 
umn density of CO and CO2 icesQ toward a sample of 
stars located behind the Taurus molecular cloud. They 

^ It is predicted that oxidation reactions involving the CO 
molecules depleted from the gas-phase can pr oduce substantial 
amounts of CO2 in the surface of dust grains l l 'Helens fc Hageiil 
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Figure 8. Pixel-by-pixel comparison between Ay and Af(CO) in the B213— L1457 filament (left) and the entire Taurus molecular cloud 
without this filament (right). 




grains is given by 



iV(CO)^ 



total 
ice 



iV(CO)icc + Af(C02)i, 



(24) 



1000 



Thus, for a given Ay the total CO column density is 
given by 



100 i:! 



N{CO) 



total 



7V(C0)gas-phasc + iV(CO) 



total 



(25) 
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Figure 9. The same as Figure |6] but including the estimated 
column density of CO and CO2 ices. For comparison we show 
the relation between visual extinction and N(CO) derived from 
observations of rare isotopic species by Frcrkine ct al. (1982) (see 
Appendix[Cll which also include the contribution for CO and CO2 
ices. 

find that the column densities are related to the visual 
extinction as 



A^(CO)i, 
10" [cm- 



and 



iV(C02)i, 



0.4(ylv-6.7), ^v >6.7mag, (22) 



0.252(^v-4.0), v4v >4.0mag. (23) 



10i7[cm-2] 
We assume that the total column of CO frozen onto dust 

[T982: Ruffic fc Herbst'l 2001l : lRos"er et al.1120011 '). Since the timescale 
of these his reactions are short compared with the cloud's lifetime, 
we need to include CO2 in order to account for the amount of CO 
frozen into dust grains along the line-of-sight. 



We can combine our determination of the column den- 
sity of gas-phase CO with that of CO ices to plot the 
total 7V(C0) as a function of Ay. The result is shown 
in Figure M The visual extinction and iV(CO)*°'^^ are 
linearly correlated over the entire range covered by our 
data, extending up to Ay — 23 mag. This result confirms 
that depletion is the origin of the deficit of gas-phase CO 
seen in Figure [S] 

In Figure [TO] we show the ratio of 7V(C0)*°''^' to 
iV(CO)gas-phase as & function of ^v, for ^v greater than 
10. The drop in the relative abundance of gas-phase CO 
from our observations is at most a factor of ~2. This is in 
agreement with previous determinations o f the depletion 
along the line of sight in molecular clouds (jKramer et al.l 
IT99alChia7eranil995l) . 

3.1.2. CO Depletion Age 

In this Section we estimate the CO depletion age (i.e. 
the time needed for CO molecules to deplete onto dust 
grains to the observed levels) in dense regions in the Tau- 
rus Molecular Cloud. We selected a sample of 13 cores 
that have peak visual extinction larger than 10 mag and 
that ^v at the edges drops below ~0.9mag (3 times the 
uncertainty in the determination of ^v)- The cores are 
located in the L1495 and B18-L1536 regions (Figure [7]). 
Unfortunately, we were not able to identify individual 
cores in Heiles's Cloud 2 due to blending. 

We first determi ne the H2 volume den sity structure of 
our selected cores. iDapp fc Basul (|2009D proposed using 
the ,Kina (,1962.) density profile, 
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Figure 10. Ratio of Af(CO)tntai to A''(CO)gas-phase plotted as 
a function of Ay for the liigli extinction portion of the Taurus 
molecular cloud. The line represents our fit to the data. 



^{r) = 



'.caVir^ 







-a') 



r <R 
r>R, 



(26) 



which is characterized by the central volume density ric, 
a truncation radius R, and by a central region of size a 
with approximately constant density. 

The column density N{x) at an offset from the core 
center x can be derived by integrating the volume density 
along a line of sight through the sphere. Defining Nc = 
2anc arctan(c) and c = R/a, the column density can be 
written 



N{x) 




2")/ arctan(c) 



1 + {x/a) 



(27) 



This column density profile can be fitted to the data. 
The three parameters to fit are (1) the outer radius R, 
(2) the central column density Nc (which in our case is 
y4v,c), and (3) the size of the uniform density region a. 

We obtain a column density profile for each core by 
fitting an elliptical Gaussian to the data to obtain its 
central coordinates, position angle, and major and mi- 
nor axes. With this information we average the data in 
concentric elliptical bins. Typical column density pro- 
files and fits to the data are shown in Figure [TTJ We 
give the derived parameters of the 13 cores we have 
analyzed in Table 2. We convert the visual extinction 
at the core center Ay.c to H2 column density assuming 
7V(H2)/Av = 9.4 X lO^^cm-^mag-^ We use then the 
definition of column density at the core center (see above) 
to determine the central volume density nc (H2) from the 
fitted parameters. 

With the H2 volume density structure, we can derive 
the CO depletion age for each core. The time needed for 
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Figure 11. Typical radial distributions of the visual extinction 
in the selected sample of cores. The solid lines represent the cor- 
responding fit. 
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Table 2 

Core Parameters 



Core ID a(J2000) (5(J2000) 



fmasl 



[pc] 



Radius 

[pc] 



nc(H2) 
[iC'cm-S] 



Mass 

[M©] 



Depletion Age 

[10^ years] 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 



04:13: 
04:17: 
04:18: 
04:18: 
04:18: 
04:19: 
04:21: 
04:23: 
04:26: 
04:29: 
04:32: 
04:33: 
04:35: 



51.63 
13.52 
05.13 
27.84 
45.66 
14.99 
08.46 
33.84 
39.29 
20.71 
09.32 
16.62 
34.29 



28:13:18.6 

28:20:03.8 

27:34:01.6 

28:27:16.3 

25:18:0.4 

27:14:36.4 

27:02:03.2 

25:03:01.6 

24:37:07.9 

24:32:35.6 

24:28:39.0 

22:42:59.6 

24:06:18.2 



22.4±0.5 
10.7±0.3 
12.3±1.2 
24.2±0.4 
9.4±0.2 
14.3±0.6 
15.2±0.3 
14.4±0.3 
15.6±0.5 
17.2±0.4 
16.0±0.5 
12.2±0.5 
12.5±0.3 



0.10±0.004 
0.19±0.021 
0.16±0.054 
0.13±0.005 
0.09±0.005 
0.12±0.010 
0.08±0.003 
0.11±0.004 
0.09±0.006 
0.13±0.006 
0.09±0.006 
0.08±0.007 
0.11±0.006 



2.01±0.30 
0.54±0.12 
0.32±0.18 
1.27±0.08 
2.00±0.93 
0.89±0.15 
1.12±0.08 
0.93±0.07 
1.48±0.23 
2.48±0.39 
3.50±2.41 
1.66±0.64 
1.64±0.35 



2.2±0.11 
0.7±0.09 
1.0±0.40 
1.9±0.07 
1.1±0.07 
1.3±0.12 
1.9±0.07 
1.4±0.06 
1.7±0.12 
1.3±0.07 
1.7±0.13 
1.5±0.14 
1.2±0.07 



307±102 
56±43 
29±63 
258±52 
110±81 
93±47 
90±19 
94±23 
143±58 
371±127 
347±347 
110±82 
145±65 



6.3±0.3 
3.4±1.5 
1.3±3.1 
3.8±0.2 
10.9±0.8 
3.1±0.6 
2.9±0.2 
5.1±0.3 
2.3±0.3 
4.5±0.3 
3.2±0.4 
6.3±0.7 
2.0±0.4 



CO molecules to deple te to a sp ecified degree onto dust 
grains is given by (e.g. lBergin fc TafaUa 2007) . 



^depiction 



5 X 109\ /n(H 



yr 



cm 



-3 



-1 
I ln(no/?igas), (28) 



where uq is the total gas-phase density of CO before 
depletion started and ngas the gas-phase CO density at 
time ^riRD ietinn- Here we assum ed a sticking cocfficienl|3 
of unity (jBisschop et al.|[2b06[ ) and that at the H2 vol- 
ume densities of interest adsorption mechanisms domi- 
nate over those of desorption (we therefore assume that 
the desorption rate is zero). 

To estimate no/ngas we assume that CO depletion oc- 
curs only in the flat density region of a core, as for larger 
radii the volume density drops rapidly. Then the total 
column density of CO (gas-phase-|-ices) in this region is 
given by N{CO)^''^ ~ 2anc(l.l x 10"''). The gas-phase 
CO column density in the flat density region of a core 
is given by iV«-*_ph_(CO) = 7V(C0)«-* - N{CO)^^-\ 

where A^(CO)i°*''' can be derived from Equation ([25|l . 
Assuming that the decrease in the [C0]/[H2] relative 
abundance in the flat region is fast and st ays constant 
toward the center of the core (models from iTafalla et alj 
p002.) suggest an exponential decrease), then no /"•gas — 
^(QQ)flat/^flat_^^^^^(CO). The derived CO depletion 
ages are listed in Table 2. Note that the fitted cores 
might not be fully resolved at the resolution of our Ay 
map (200"or 0.14pc at a distance of 140pc). Although 
rio/rigas is not very sensitive to resolution, due to mass 
conservation, we might be underestimating the density 
at the core center. Therefore, our estimates of the CO 
depletion age might be considered as upper limits. 

In Figure [T2] we show the central density and the cor- 
responding depletion age of the fitted cores as a function 
of ^v The central volume density is well correlated 
with Ay but varies only over a small range: its mean 
value and standard deviation are (1.4 ± 0.4)xl0*cm~^. 
Still, the moderate increase of n(H2) with Ay com- 
pensates for the increase of iV (CO)totai/A^(CO)gas-phasc 
with Ay to produce an almost constant depletion age. 
The mean value and standard deviation of idcpiction are 

* The sticking coefficient is de fined as how often a spe cies will 
remain on the grain upon impact IjBergin fc TafalTall2007l ). 
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Figure 12. {upper panel) The central H2 volume density as a 
function of the peak Av for a sample of 13 cores in the Taurus 
molecular cloud. The line represent a fit to the data, {lower panel) 
CO depletion age as a function of Av for the sample of cores. 
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(4.2 ±2.4) X 10^ years. This suggests that dense cores at- 
tained their current central densities at a similar moment 
in the history of the Taurus molecular cloud. 

3.2. Low N(^^CO) column densities 

In the following we compare the lowest values of the 
CO column density in our Taurus survey with the visual 
extinction derived from 2MASS stellar colors. In Fig- 
ure [13] we show a comparison between A^(CO) and Ay 
for values lower than 5 magnitudes of visual extinction. 
The figure includes CO column densities for pixels lo- 
cated in Mask 1 and 2. We do not include pixels in Mask 
because its single value does not trace variations with 
Ay. Instead, we include a horizontal line indicating the 
derived average CO column density in this Mask region. 
We show a straight line (blue) that indicates N{CO) ex- 
pected f rom a abundance ratio [^^CO]/[H2]=l.lxlO"'' 
(Section 13.11) . The points indicate the average Ay in a 
A^(CO) bin. We present a fit to this relation in Figure [HI 

The data are better described by a varying [^^C0]/[Il2] 
abundance ratio than a fixed one. This might be caused 
by photodissociation and fractionation of CO which can 
produce strong variations in the CO abundances between 
UV-exposed and shielded regions (van Dishocck & Blac3 
[19881: [Visser et al.[[2009[ ). To test this possibility we in- 
clude in the figure severa l models of these effects pro- 
vided by Ruud Visser f see [Visser et al.[[2009l for details). 
They show the relation between Ay and A^(H2) for dif- 
ferent values of the FUV radiation field starting from 
X = 1.0 to 0.1 (in u nits of the me an interstellar radia- 
tion field derived by [Draind [1978D . All models have a 
kinetic temperature of 15K and a total H volume den- 
sity of 800 cm~^ which corresponds to n(H2) ~ 395 cm~'^ 
assuming n(IIi)= 10cm~^. (This value of n(H2) is close 
to the average in Mask 1 of 375 cm~^.) The observed re- 
lation between Ay and A^(CO) cannot be reproduced by 
a model with a single value of x- This suggests that the 
gas have a range of physical conditions. Considering the 
average value of Ay within each bin covering a range in 
A^(CO) of 0.25 dex, we see that for an increasing value 
of the visual extinction, the FUV radiation field is more 
and more attenuated so that we have a value of A^(CO) 
that is predicted by a model with reduced x- 

We also include in F i gure [13[ the fit to the observa- 
tions from [Sheffer et al.l ([200^ toward diffuse molecular 
Galactic lines-of-sight for log(7V(H2)) > 20.4. The fit 
seems to agree with the portion our data points that 
agree fairly well wi th the model having x = l-O- Since 
[Sheffer et al.[ ([2008[ ) observed diffuse lines-of-sight, this 
suggests that a large fraction of the material in the Tau- 
rus molecular cloud is shielded against the effect of the 
FUV illumination. T his is supporte d by infrared obser- 
vations in Taurus by iFlagev et al.l ([2009) that suggest 
that the strength of the FUV radiation field is between 
X = 0.3 and 0.8 . 

[Sheffer et all ([2008[ . see also [Federman eFall [198O0 
showed empirical and theoretical evidence that the scat- 
ter in the Ay — N{CO) relation is due to variations of 
the ratio between the total H volume density (n^'^' = 
riH + 2nH2 ) ^^^ the strength of the FUV radiation field. 
The larger the volume density or the weaker the strength 
of the FUV field the larger the abundance of CO relative 
to II2. Note that the scatter in the observations from 



[Sheffer et al.[ ([2008[ ) is much smaller than that shown in 
Figure [T3| This indicates that we are tracing a wider 
range of physical conditions of the gas. The exc i tation 
temperature of the gas observed by iSheffer et al.l (|2008f ) 
does not show a large variation from Tcx=5 K while we 
observe values between 4 and 15 K. 

In Figure [T3| we see that some regions can have large 
[^^C0]/[Il2] abundance ratios but still have very small 
column densities {Ay = 0.1 — 0.5 mag). This can be 
understood in terms of a medium which is made of an 
ensemble of spatially unresolved dense cluinps embedded 
in a low density interclump medium ([BenschI 2006). In 
this scenario, the contribution to the total column den- 
sity from dense clumps dominates over that from the 
tenuous inter-clump medium. Therefore the total col- 
umn density is proportional to the number of clumps 
along the line-of-sight. A low number of clumps along a 
line-of-sight would give low column densities while in the 
interior of these dense clumps CO is well shielded against 
FUV photons and therefore it can reach the asymptotic 
value of the [C0]/[H2] ratio characteristic of dark clouds. 

4. DISCUSSION 
4.1. The mass of the Taurus Molecular Cloud 

In this section we estimate the mass of the Taurus 
Molecular Cloud using the A^(CO) and Ay maps. The 
masses derived for Mask 0, 1, and 2 are listed in Table 3. 
To derive the II2 mass from iV(CO) we need to apply an 
appropriate [C0]/[H2] relative abundance for each mask. 
The simplest case is Mask 2 where we used the asymp- 
totic i^CO abundance of 1.1x10"'' (see Section [XTjl . We 
corrected for saturation including temperature gradients 
and for depletion in the mass calculation from 7V(C0). 
These corrections amount to ~319M0 (4Mq from the 
saturation correction and 315 Mq from the addition of 
the column density of CO-ices). For Mask 1 and 0, we 
use the fit to the relation between iV(Il2) and A^(CO) 
shown in Figure [T4[ As we can see in Table 3, the masses 
derived from Ay and iV(CO) are very similar. This con- 
firms that N{CO) is a good tracer of the bulk of the 
molecular gas mass if variations of the [C0]/[Il2] abun- 
dance ratio are considered. 

Most of the mass derived from Ay in Taurus is in Mask 
2 (~49%). But a significant fraction of the total mass lies 
in Mask 1 (-28%) and Mask (-23%). This imphes that 
mass estimates that only consider regions where ^'^CO is 
detected underestimate the total mass of the molecular 
gas by a factor of —2. 

We also estimate t he masses of high-colum n density 
regions considered by [Goldsmith et al.' f2008') that were 
previously defined bylOnishi et al. (1996). In Table 4 we 
list the masses derived from the visual extinction as well 
as from N{CO). Again, both methods give very similar 
masses. These regions together represent 43% of the total 
mass in our map of Taurus, 32% of the area, and 46% 
of the ^^CO luminosity. This suggest that the mass and 
^^CO luminosity are uniformly spread over the area of 
our Taurus map. 

A commonly used method to derive the mass of molec- 
ular clouds when only ^^CO is available is the use 
of the empirically derived CO-to-H2 conversion factor 
(Xco = iV(H2)//co - MujLco )• Observations of 7- 
rays indicate that this factor is 1.74xlO^'^cm~^(K km 
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Figure 13. Comparison between tiie visual extinction derived from 2MASS stellar colors and the ^ CO column density derived from 
^•^CO and ^^CO observations in Taur us for Ay < 5 mag . The blue line represents the ^^CO column density derived from A\r assuming 
^(H2)/Av = 9-4 X lO^^cm"^ mag~^ IBohlin et al.in~97a) and a [CO]/[H2] abundance ratio of 1.1 X 10""*. The gray scale represents the 
number of pixels of a given value in the parameter space and is logarithmic in the number of pixels. The red contours are 2,10,100, and 
1000 pixels. The black lines represent several models of selective CO photodissociation and fractionation provided by Ruud Visser (see 
text). The light blue line represents the fit to the observations from Sheffer et al. (2008) toward diffuse molecular Galactic lines-of-sight 
for log(Af(H)2) > 20.4. The horizontal line represents the average N{CO) derived in Mask 0. Each pixel has a size of 100" or 0.07 pc at a 
distance of 140 pc. 

s-1 pc-2 W^M(M„)=3.7ico (K Km s'^ pc^) in our 
Galaxy (jGrenier et a l. 2005). To estimate Xqo in Mask 
2, 1, and we calculate the ^^CO luminosity (Leo) in 
these regions and compare them with the mass derived 
from Ay. We also calculate Xco from the average ratio 
of iV(H2) (derived from Ay ) to the CO integrated inten- 
sity Ico for all pixels in Maskl and 2. For MaskO, we 
used the ratio of the average A^(H2) (derived from Ay 
) to the average CO integrated intensity obtaining after 
combining all pixels in this mask region. The resulting 
values are shown in Table 3. The table shows that the 
difference in Xco between Mask 2 and Mask 1 is small 
considering that the [C0]/[H2] relative abundance be- 
tween these regions can differ by up to two orders of 
magnitude. The derived values are close to that found in 
our Galaxy using 7-ray observations. For MaskO, how- 
ever, Xco is about an order of magnitude larger than in 
Maskl and 2. 

Finally we derive the surface density of Taurus by com- 
paring the total H2 mass derived from Ay (15015 M©) 
and the total area of the cloud (388 pc'^). Again, we 
assumed that in Mask the CO-emitting region occu- 
pies 26% of the area. The resulting surface density is 
~39MqPC~^ which is very similar to the median value 
of 42 Mq pc~^ derived from a large sample of galactic 




Figure 14. The average A'^(H2) and Ay as a function of N{CO) 
in Maskl. N(R2) is estimated from Ay assuming N{}l2)/Ay = 



9.4 X lO^Oc 
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Figure 15. Probability density function of the visual extinc- 
tion in the Taurus molecular cloud. The solid line corresponds 
to a Gaussian fit to the distribution of the natural logarithms of 
Av/{Av). This fit considers only visual extinctions that are lower 
than 4.4 mag, as the distribution deviates clearly from a Gaussian 
for larger visual extinctions, (see text). 

molecular clouds by iHever et al.l ()2009D . 

4.2. Column density probability density function 

Numerical simulations have shown that the probabil- 
ity density function (PDF) of volume densities in molec- 
ular clouds can be fitted by a log-normal distribution. 
This distribution is found in simulations with or with- 
out magnetic fields when self-gravity is not important 
JOstriker et al. 200l HNordlund fc PadoanlUQQalLTet al.l 
l2004t iKlessenii2000l) . A log- normal distribution arises as 
the gas is subject to a succession of independent compres- 
sions or rarefactions that produce multiplicative varia- 
tions of the volume density ( Passo t & Vaz q uez- Semadenil 
119981: iVazauez-Semadeni fc G~c ia 2001). This effect 
is therefore additive for the logarithm of the vol- 
ume density. A log-normal function can also de- 
scribe the distribution of column densities in a molec- 
ular cloud if compressions or r a.refactions along the 
line of sight are independent ( Ostriker et al.l 120011 : 
IVazauez-Semadeni fc Garcial 120011 )^ Note that log- 



normal distributions are not an exclusive result of su- 
personic turbulence as they are also seen in simulations 
with the presence of self-gravity and /or strong ma gnetic 
fields but without strong turbulence (jTassis et al.l[2010( ). 

Deviations from a log-normal in the form of 
tails at high or low densities are expected if the 
equation of state deviates from being isotherma l 
(jPassot fc Vazauez-Semadenil 119981: iScalo et al.l I1998D . 
This, however, also occurs in simulations with an isother- 
mal equation of st ate due to the effects of self-gravity 
(jTassiset al.llMol ). 

In Figure [15] we show the histogram of the natural log- 
arithm of Ay in the Taurus molecular cloud normalized 
by its mean value (1.9 mag). Defining x = N/{N), where 
N is the column density (either Ay or N{Il2) ), we fit a 
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Figure 16. Probability density function of the visual extinction 
for Maskl {upper panel) and Mask 2 (lower panel) in the Taurus 
molecular cloud. The solid line corresponds to a Gaussian fit to 
the distribution of the natural logarithm of Ay /{Ay). The fit 
for Maskl considers only visual extinctions that are larger than 
0.24 mag, while the fit for Mask 2 includes only visual extinctions 
that are less than 4.4 mag (see text). 



function of the form 
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A^pixois exp 
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(29) 



where ^ and a^ are the mean and variance of ln(a;). The 
mean of the logarithm of the normalized column den- 
sity is related to the dispersion cr by /i = — ct^/2. In all 



Gaussian fits, we consider y/N counting errors in each 
bin. 

The distribution of column densities derived from the 
visual extinction shows tails at large and small Ay. The 
large- Ay tail starts to be noticeable at visual extinc- 
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Table 3 

Properties of Different Mask Regions in Taurus 



Region # of Pixels" 



Mass from 
i^CO and i^co 

[Mq] 



Mass from Area 



[CO]/[H2 



[Me] 



[PC 



Leo ^£o = ^(H2)//co ^So = ^/^co 

[K km s-i pc2] [cm-2/(K km s-^)] [Mq/{K km s"! pc^) 



Mask 


52338 


3267 


3454 


63'' 


1.2x10-'' 


130 


1.2x10^1 


Mask 1 


40101 


3942 


4237 


185 


variable 


1369 


1.6x10^0 


Mask 2 


30410 


7964 


7412 


140 


1.1x10-4 


1746 


2.0x1020 


Total 


122849 


15073 


15103 


388 




3245 


2.3xl020 



26 
3.1 

4.2 
4.6 



^ At the 200^^ resolution of the A\r map. 

Calculated from the average ratio of A'^(H2), derived from Av, to CO integrated intensity for each pixel. 
^ Total mass per unit of CO luminosity. 

Effective area of CO emission based in the discussion about Mask in Seetion l2.3l 



Table 4 

Mass of Different High Column Density Regions in Taurus 



Region 


# of Pixels 


Mass from 
13C0 and i2cO 


Mass from 


Area 


^co 






[Mr,] 


[Mq] 


[pc^] 


[K km s-i pc2] 


L1495 


7523 


1836 


1545 


35 


461 


B213 


2880 


723 


640 


13 


155 


L1521 


4026 


1084 


1013 


19 


236 


HCL2 


3633 


1303 


1333 


17 


221 


L1498 


1050 


213 


170 


5 


39 


L1506 


1478 


262 


278 


7 


68 


B18 


3097 


828 


854 


14 


195 


L1536 


3230 


474 


579 


15 


134 


Total 


26917 


6723 


6412 


125 


1509 



Nu [102lcm-2] 
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Figure 17. Probability density function of the H2 column den- 
sity derived from N(CO) in Mask 2 with an angular resolution of 
47" (0.03 pc at the distance of Taurus, 140 pc). The solid line corre- 
sponds to a Gaussian fit to the distribution of the natural logarithm 
of Af(H2)/(A'^{H2)). The fit considers H2 column densities that are 
lower than 4x10^^ cm— ^ (or ~4mag). 

tions larger than ~4.4mag. The low- Ay tail starts to be 
noticeable at visual extinctions smaller than ~ 0.26 mag, 
which is similar to the uncertainty in the determination 
of visual extinction (0.29 mag), and therefore it is not 



possible to determine whether it has a physical origin 
or it is an effect of noise. The distribution is well fit- 
ted by a log-normal for Ay smaller than 4.4 mag. We 
searched in our extinction map for isolated regions with 
peak Ay > 4.4 mag. We find 57 regions that satisfy this 
requirement. For each region, we counted the number 
of pixels that have Ay > 4.4 mag and from that calcu- 
lated their area, A. We then determined their size using 
L — 2y/(A/Tr). The average value for all such regions 
is 0.41 pc. This value is similar to the Jeans length, 
which for Tkin=10K and n(H2) = 10^ cm-'^ is about 
0.4 pc. This agreement suggests that the high- Ay tail 
might be a result of self- gravity acting in dense regions. 
iKainulainen et al] ()2009[ ) studied the column density dis- 
tribution of 23 molecular cloud complexes (including the 
Taurus molecular cloud) finding tails at both large and 
s mall visual extincti o ns. 

IKainulainen et all ()2009( ) found that high- Ay tails are 
only present in active star-forming molecular clouds 
while quiescent clouds are well fitted by a log-normal. 
We test whether this result applies to regions within 
Taurus in Figure [16] where we show the visual extinc- 
tion PDF for the Mask 1 and 2 regions. Mask 1 includes 
lines-of-sights that are likely of lower volume density 
than regions in Mask 2, and in which there is little star 
formation. This is illustrated in Figure [T] where we show 
the distribution of the Mask regions defined in our map 
ov erlaid by the co mpilation of stellar members of Taurus 
bv lLuhman et all (|^006). Most of the embedded sources 
in Taurus are located in Mask 2. Note that the normal- 
ization of Ay is different in the two mask regions. The 
average value of Ay in Maskl is 0.32 mag and in Mask 2 
is 2.1 mag. In Maskl we see a tail for low- Ay starting 



The relation between gas and dust in the Taurus Molecular Cloud 



17 



at about 0.2 mag. Again, this visual extinction is close to 
the uncertainty in the determination of Ay. For larger 
visual extinctions the PDF appears to be well fitted by a 
log-normal distribution. In case of the visual extinction 
PDF in Mask 2, we again see the tail at large Ay starting 
at about 4.4 mag. For lower values of Ay the distribution 
is well represented by a log-normal. 

We can use our CO map of Taurus at its original reso- 
lution (47" which corresponds to 0.03 pc at the distance 
of Taurus, 140 pc) to study the column density PDF at 
higher resolution than the 200" Ay map fFigure [T7| . We 
estimate A^(H2) from our CO column density map in 
Mask 2 by applying a constant [C0]/[H2] abundance ra- 
tio of 1.1x10^^ (Section 13. ip . The average H2 column 
density in Mask 2 is 3 x lO^^cm"^. We do not con- 
sider Maskl because of the large scatter found in the 
[C0]/[H2] abundance ratio (Section 13. 2p . In the figure 
we see that the distribution is not well fitted by a log- 
normal. As for Ay, the PDF also shows a tail for large 
column densities that starts to be noticeable at about 
4 X 10^^ cm^^ (or Ay ~ 4 mag). Therefore, the high- 
column density excess seems to be independent of the 
spatial scale at which column densities are sampled. We 
repeated the procedure described above to search for iso- 
lated cores in our map with A^(H2) > 4 x 10^^ cm~^ and 
obtained an average size for cores of 0.5 pc, which is con- 
sistent to that obtained in our Ay map . Note that at this 
resolution we are not able to account for effects of tem- 
perature gradients and of CO depletion along the line of 
sight, as this requires knowledge of Ay at the same reso- 
lution. We therefore underestimate the number of pixels 
in the H2 column density PDF for A^(H2) > 1 x 10^^ cm~^ 
while we overestimate them for A^(H2) < 1 x 10^^ cm~^. 
But the number of pixels (^-^7000) affected by those ef- 
fects represent only 9% of the number of pixels (^81000) 
that are in excess relative to the log-normal fit between 
3 X 10^"'^ and 1 x lO^^cm"^, and therefore the presence 
of a tail at large- A^(Il2) is not affected. Note that this 
also affected our ability to identify isolated regions in 
the A^(H2) map. We were able to identify only 40 cores 
compared with the 57 found in the ^v niap. 

In summary, we find that the distribution of column 
densities in Taurus can be fitted by a log-normal distri- 
bution but shows tails at low and high-column densities. 
The tail at low-column density may be due to noise and 
thus needs to be confirmed with more sensitive maps. 
We find that the tail at large column densities is only 
present in the region where most of the star formation 
is taking place in Taurus (Mask 2) and is absent in more 
quiescent regions (Maskl). The sa me trend has been 
found in a larger sample of clouds by iKainulainen et al.l 
()2009() . Here we suggest that the distinction between 
star-forming and non star-forming regions can be found 
even within a single molecular cloud complex. The pres- 
ence of tails in the PDF in Taurus appears to be inde- 
pendent of angular resolution and is noticeable for length 
scales smaller than 0.41 pc. 

5. CONCLUSIONS 

In this paper we have compared column densities de- 
rived from the large scale ^^CO and ^' ^CO maps of the 
Taurus molec ular cloud presented by iNaravanan et al.l 
120081 (see also lGoldsmith et al.|[2008l ) with a dust extinc- 
tion map of the same region. This work can be summa- 



rized as follows, 

• We have improved the derivation of the CO 
column density compared to that derived by 
[Goldsmith et al.l (HOOs) by using an updated value 
of the spontaneous decay rate and using exact nu- 
merical rather than approximate analytical calcu- 
lation of the partition function. We also have used 
data that has been corrected for err or beam pick- 
up us ing the method presented by iBensch et al.l 
(I2OOI . 

• We find that in the Taurus molecular cloud the col- 
umn density and visual extinction are linearly cor- 
related for Ay up to 10 mag in the region associated 
with the B213-L1495 filament. In the rest of Tau- 
rus, this linear relation is fiattcned for Ay > 4 mag. 
A linear fit to data points for Ay < 10 mag and 
A^(CO) > lO-'^'^cm"^ results in an abundance of 
CO relative to H2 equal to 1.1x10"^. 

• For visual extinctions larger than ~4mag the CO 
column density is affected by saturation effects and 
freezeout of CO molecules onto dust grains. We 
find that the former effect is enhanced due to the 
presence of edge-to-center temperature gradients 
in molecular clouds. We used the RATRAN ra- 
diative transfer code to derive a correction for this 
effect. 

• We combined the column density of CO in ice form 
derived from observatio ns towards embedde d and 
field stars in Taurus bv lWhittet et al.l (|2007[ ) with 
the saturation-corrected gas-phase 7V(C0) to de- 
rive the total CO column density (gas-phase-|-ices). 
This quantity is linearly correlated with Ay up to 
the maximum extinction in our data ^23 mag. 

• We find that the gas-phase CO column density is 
reduced by up to a factor of ~2 in high-extinction 
regions due to depletion in the Taurus molecular 
cloud. 

• We fit an analytical column density profile to 
13 cores in Taurus. The mean value and stan- 
dard deviation of the central volume density are 
(1.4 ± 0.4)xl0^cm~'^. We use the derived volume 
density profile and the amount of depletion ob- 
served in each core to derive an upper limit to the 
CO depletion age with a mean value and standard 
deviation of (4.2 ± 2.4) x 10^ years. We find lit- 
tle variation of this age among the different regions 
within Taurus. 

• For visual extinctions lower than 3 mag we find that 
iV(CO) is reduced by up to two orders of magni- 
tude due to the competition between CO formation 
and destruction processes. There is a large scatter 
in the Ay -iV(CO) relation that is suggestive of dif- 
ferent FUV radiation fields characterizing the gas 
along different lines-of-sight. 

• The mass of the Taurus molecular cloud is about 
1.5xlO'*Mo. Of this, ^49% is contained in pixels 
where both ^^CO and ^^CO are detected (Mask 2), 
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-28% where ^^CO is detected but ^^CO is not 
(Mask 1), and -23% where neither ^^CO nor i^CO 
are detected (MaskO). 

• We find that the masses derived from CO and Ay 
are in good agreement. For Mask 2 and MaskO we 
used a [C0]/[H2] relative abundance of 1.1 xlO^'' 
and 1.2x10^^, respectively. For Maskl, we used a 
variable [C0]/[H2] relative abundance taken from a 
fit to the average relation between Ay and iV(CO) 
in this region, with -6.7 < log([CO]/[H2]) < -3.9. 

• We also compared the mass derived from Ay with 
the ^^CO J = 1 — >■ luminosity for the regions de- 
rived above. For Mask 1 and 2 these two quantities 
are related with a CO-to-H2 conversion factor of 
about 2.1xl02°cm-2 (K km s'^ )-^. The derived 
CO-to-H2 conversion factor is in agreement with 
that found in our Galaxy using 7-ray observations. 
In MaskO, however, we find a larger the conversion 
factor of 1.2xl02icm-2 (K km s^^ y^. 

• We studied the distribution of column densities in 
Taurus. We find that the distribution resembles 
a log-normal but shows tails at large and low 
column densities. The length scale at which the 
high-column density tail starts to be noticeable 
is about 0.4 pc, which is similar to the Jeans 
length for a T=10K and tihs = lO'^cm"'^ gas, 
suggesting that self-gravity is responsible for its 
presence. The high-column density tail is only 



present in regions associated with star formation, 
while the more quiescent positions in Taurus do 
not show this feature. This tail is independent of 
the resolution of the observations. 
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APPENDIX 
ERROR BEAM CORRECTION 

The FCRAO 14m telescope is sensitive not only to emission that couples to the main beam (with efiiciency ?7mb = 0.45 
at 115 GHz and ?7mb = 0.48 at 110 GHz) but also to emission distributed on scales comparable to the error beam 
(30'). For emission extended over such large-scale, the coupling factor (including the main beam contribution) is the 
forward spillover and scattering efficiency, Tyfss = 0.7, at both frequencies. The error beam pickup, also known as "stray 
radiation" , can complicate the accurate calibration of the measured intensities: a straightforward scaling of the data 
by l/yymb can significantly overestimate the true intensity in regions where emission is present on large angular scales. 
Given the wide range of angular sizes of the structures in Taurus, it is clear in general that neither ?7nib nor rjfss will 
give optimum results. 

To accurately scale the FCRAO data, it is essential to remove the error beam component before scaling the intensities 
to the main beam scal e. Methods for correcting millimeter-wave data for error beam pickup have been discussed 
bv iBensch et al.l (|200lD . who introduce the "corrected main beam temperature scale" (T,„b,c) with which optimum 
calibration accuracy is achieved by scaling the data by l/r/,„b after removal of radiation detec ted by the e rror beam. 

To remove the error beam component we use the second of the methods described in Bensch et al.l (|2001l ). The error 
beam component is removed in Fourier space directly from the FCRAO data, and the intensities are converted to the 
Tlnb.c scale, by the following method: 



• The Fourier transform of the antenna temperature is taken: TJ = FT{Tp^) 

• The following correction is applied to each velocity slice in the cube: 



Tinb.c = Tl{i]nib + rich exp(- 



-n^e! 



41n(2) 



-'mhl 



))- 



where O^h is the FWHM of the error beam, 6'nib is the FWHM of the main beam, 770b = Viss — ?ymb, and kx, ky 
are the wavenumbers along the x, y (RA, deck) directions respectively. 

• The inverse Fourier transform is performed, with only the real part of the result being retained: 
Tmh,c — Re (/i^T(r,nb,c))- Thc imaginary part is consistent with round-off errors. 

At low spatial frequencies, the correction factor is — l/'7fss while at high spatial frequencies, the correction factor is 
' l/?7mb- The effective correction factor at any point is determined therefore by the spatial structure of the emission 
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Figure 18. fleft) Color (H — Kg) versus (J — H) for stars observed in the control field, (right) Intrinsic J — H and H — Kg colors of 
Main Sequence, Giant, and Supergiant stars (taken from IKoornneeJI 1983) . These stars are indicated by circles in the left panel. 



in the vicinity of that point. More detailed information and quantitative analysis of the above procedure can be 
found in Brunt et al (2010, in prep) and Mottram & Brunt (2010, in prep). For typical applic ations, a naive sca ling 
by l/?7mb overestimates the true intensities, as inferred from comparison to CfA survey data (jDame et al.ll200ll ). by 
around 25-30%. For reference, an overestimation of ~ 50% would be applicable if rjiss were appropriate everywhere. 
The spatially variable correction factor afforded by the method used here therefore offers a higher fidelity calibration 
of the data. 

THE EXTINCTION MAP 

We have used the 2MASS point s ource catalog to create an near -infrared extinct ion map of Taurus. This was done 
using an implementation of NICER (jLombardi fc Alvesl[2001[ ) from lChapmanI ()2007[) . The 2MASS catalog we used has 
1039735 (--1 million) stars over an area bet ween RA= 04:03:51.6 and 05:05:56.6 and decl.=+19:24:14.4 and -f 30:50:24 
(J2000). We use the compilation by Luhm an et al.l |2006) to remove 156 stars that are known to be members of 
Taurus. The map generated has an angular resolution of 200" and is Nyquist sampled with a pixel spacing of 100", 
corresponding to 0.07 pc at a distance of 140 pc. The resolution of the map was determined by that of the Hi map 
used to correct the data for the contribution of H i to the total extinction (see below) . The final extinction map is 
shown in Figure 2] 

We constructed extinction maps in nearby regions around Taurus with the goal of finding a field that does not show 
significant extinction, so it can be used as a control field to estimate the intrinsic (J — H) and {H — Kg) stellar colors. 
We selected a region corresponding to a 2° x 2° box centered at RA=03:50:44.7 and decl.=-t-27:46:54.1 (J2000). The 
mean (± weighted standard deviation) values for stars in this box are 0.454±0.157 mag for (J — H) and 0.114±0.074 
mag for {H — Kg). We also computed the covariance matrix for the {J — H) and {H — Kg) colors. The on-axis elements 
of this matrix are (Tj_^ and crj^^x i ^he dispersions of the (J — H) and {H — Kg) colors in the control field, while 
the two off-axis elements are identical to each other, with a value of 0.006. In Figure [18] we show the color {H — Kg) 
versus (J — H) of stars in the control field. We also show the intrinsic color of Main Sequence, Giant, and Supergiant 
stars. Apart from the scatter due to photometric errors, there is large scatter in the intrinsic colors due to different 
stellar types in the control field. The mean values, weighted standard deviation and off-axis covariance matrix are 
input to the NICER ro utine and with them we correct for the different sources of scatter of the intrinsic colors in the 
control field. Note that iPadoan et all ()2002| ) used an intrinsic {H — Kg) color of 0.13 mag in their extinction map of 
Taurus. The difference relative to that in our control field is thus seen to be small. 

We transformed the (J — H) and {H ~ Kg) colors to Ay using an extinction curve from iWeingartner fc Draiiii 
()2001|) with a ratio of selective to total extinction i?v=3-l. Note that this R\ is derived towards diffuse regions (Ay < 
1.5 mag). At larger volume densities the value oi Ry is expected to increase up to 4.5 at the center of dense cores 
due to grain growth by accretion and coagulation ("Whitte t et al.l[2001l ). Considering that a given line-of-sight might 
intersect both dense and diffuse regions, iWhittet et a l. (200l|) estimated an effective i?v that increases up to ~ 4.0 for 
Ay — 10 mag. For such a value of Ay we expect that for a given total hydrogen column dens ity, the visual extinction 
increases by about 20% due to enhanced scattering as the grain sizes increases ()Drainell2003f) . 

Figure [12] shows the relation between the estimated visual extinction in the Taurus molecular cloud and the formal 
error per pixel (i.e. error propagation from the error in the estimation of Ay for each star) and the number of stars 
per pixel. The errors in Ay range from ~0.2mag at low extinctions to ~1.3mag at large visual extinctions. The 
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Figure 19. The formal error per pixel (left) and the number of stars per pixel (right) as a function of the visual extinction estimated in 
the Taurus molecular cloud. 
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Figure 20. Histogram of the visual extinction in Taurus associated with Hi for the full range of 21cm velocities (light gray) and for the 
velocity range between and 20 km s~^ (dark gray). The histogram has been determined for the region where both ^^CO and ^■'CO are 
detected. 



average error is 0.29 mag while the average number of stars per pixel is 28. As expected, the number of stars per pixel 
decreases as the extinction increases. 

There is a large filamentary Hi structure, extending away from the Galactic plane, which coincides with the eastern 
part of Taurus. Based on distances of molecular clouds at the end of the filament, we assume that this filament lies 
between the Taurus background stars and the Earth. Dust in the filament will thus contribute to the total extinction 
measured. We estimate the contribution to the visual extinction from dust associated with H i using the Arecibo map 
from Marco Krco (PhD Thesis, Cornell University, in preparation). In Figure [20l we show a histogram of the visual 
extinction associated with positions in the Hi map where both ^^CO and ^'^CO are detected in our Taurus map. We 
show the extinction for the full range of velocities and for the range between to 20 km s^^ (similar to the velocity 
range at which CO emission is observed). We correct the Ay map by extinction associated with neutral hydrogen in 
the latter range (see below). The average correction is ~0.3mag. 

In order to see whether some H i velocity components are foreground to the 2MASS stars, we examine a field with 
complex Hi velocity structure. We choose a region northwest of the Taurus molecular cloud that shows small visual 
extinction (RA = 04:57:25.472 and deck = 29:07:0.81). Figure [2Tk shows a histogram of the visual extinction without 
correction, corrected for Hi over the entire velocity range, and corrected for Hi over the to 20 km s~^ range. We 
also show in Figure [2Tb the average H i spectrum in the selected field. The negative velocity components produce 
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Figure 21. (left) Histogram of the visual extinction calculated in a field with complex velocity structure (see text) with and without 
correction for extinction associated with Hi. (right) The Hi spectrum averaged over the this field. 



significant excess reddening associated with Hi that is inconsistent with the extinction determined from 2MASS stars. 
We therefore conclude that Hi components with negative velocities are background to the 2MASS stars. This confirms 
the correctness of excluding negative velocities for determining the H l-associated extinction correction for Taurus. 
The exact velocity range used is a source of uncertainty of a few tenths of a magnitude in the extinction. 

We finally note that the widespread H i emission is also present in the control field. The control field is contaminated 
by ^0.12 mag of visual extinction associated with Hi. This contribution produces a small overestimation of the 
intrinsic colors in the control field. Therefore, since we determine visual extinctions based on the difference between 
the observed stellar colors in Taurus and those averaged over the control field, we have added 0.12 mag to our final 
Ay map of Taurus. 

CORRECTION FOR TEMPERATURE GRADIENTS ALONG THE LINE-OF-SIGHT 

In order to assess the impact of core-to-edge temperature gradients in the estimation of A^(CO), we use the radiative 
transfer code RATRAN ()Hogerheiide fc van der Takll200dl) . With RATRAN we calculate ^^CO and ^^CO hue profiles 
and integrated intensities from a model cloud and use them to estimate iV(CO), following our analytic procedure 
(Section [2]) , which we then compare with that of the original model cloud. 

The model is a spherical cloud with a truncated power-law density profile. We adopt a density profile 
n(r)=ns(r/rc)~" for 0.3rc < r < Tc, and constant density, n(r) — 71,^(0.3)"" in the central portion of the cloud 
(r < 0.3rc). Here, re is the cloud radius and Us the density at the cloud surface. In our models we use a power-law 
exponent of a = 1.96 for the density profile, a density at the cloud surface of 2x10'' cm~^, and a cloud radius of 
5x10^^ cm (O.OlSpc at a distance of 140pc). The values of the density at the cl oud surface and the power-law 
exponent are taken from the fit to the extinction profile of the prestellar core B68 (lAlves et al.ll200il) . We adopt a 
FWHM line-width in the model of Au=1.33 km s""', in order to have line widths that are consistent with those in 
pixels having Ay > 10 mag f Figure [22]) . We use the ^^CO and ^^CO emission resulting from the radiative transfer 
calculations to derive, using Equation ([TT]). the CO column density (Af(CO)omission) that will be compared with that 
of the model cloud (A^(CO)modci)- We trace CO column densities between 3xl0^^cm~^ and IxlO-'^cm^^. 

We consider the case of isothermal clouds and of clouds with temperature gradients. Temperature gradie nts can 
be pr oduced, for example, when clouds are externally illuminated by the interstellar radiation field (e.g. lEvans et al.l 
1200 It) . We adopt a temperature profile given by 



T{r) - (T, - T,) 



Tr. 



(CI) 



where T^ and Tg are the temperature at the cloud center and surface, respectively. 

We run isothermal cloud models with kinetic temperatures of 8, 9, 10,12, and 15 K. In the case of clouds with 
temperature gradients we consider the same range of temperatures for the cloud surface and Tc = 8 K for the cloud 
center. We choose this value because the balance between the dominant heating and cooling mechanisms in dense and 
shielded regions, namely cosmic-ray heating and cooling by gas-grain collisions, typically results in this range of tem- 
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peratures ()Goldsmithll200lD . The selected range of temperatures match the observed range of excitation temperatures 
for Ay > 10 mag (Figure [22]) . 
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Figure 22. Pixcl-by-pixel comparison between the visual extinction (^v) and the ^■^ CO-derived excitation temperature (left 
FWHM line-width (right). 
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In the upper rows of Figure [^ we show the excitation temperature derived from the model ^^CO emission as a 
function of iV(CO)modci- The left panel corresponds to isothermal clouds and the right panel to clouds with temperature 
gradients. We also show the line-of-sight (LOS) averaged kinetic temperature as a function of -/V(CO),nodci- In clouds 
with temperature gradients, low column densities are on average warmer than larger column densities. For isothermal 
clouds Tox and the model kinetic temperature are almost identical for all values of iV(CO)inodci. In the case of clouds 
with temperature gradients we see that, although the average LOS kinetic temperature decreases for large CO column 
densities, the derived excitation temperature shows little variation with A^(CO)modci, tracing only the temperature at 
the cloud surface. This is a result of ^^CO becoming optically thick close to the cloud surface and therefore the T^x 
determined in this manner applies only to this region. In clouds with temperature gradients, using ^^CO to calculate 
the excitation temperature overestimates its value for regions with larger column densities, where most of the ^^CO 
emission is produced. 

We show the ^^CO opacity (Equation [20]) for both isothermal clouds and clouds with temperature gradients as a 
function of iV(CO)modei in the middle panels of Figure [23l We also show opacities calculated from the model cloud 
^■^CO column densities assuming LTE (tlte)- When the cloud kinetic temperature is constant, both opacities show 
good agreement for all sampled values of iV(CO)modci. In contrast, for clouds with temperature gradients, opacities 
derived from the model line emission are lower than tlte by up to a factor of ^3. The differences arise due to the 
overestimation of the excitation temperature in regions with large iV(CO), as Equation |20| assumes a constant value 
of Tex. 

The relation between A^(CO)emission and A^(CO)modoi is shown in the lower panels of Figure [131 Isothermal clouds 
show almost a one-to-one relation between these two quantities whereas clouds with temperature gradients show that 
the relation deviates from linear for large iV(CO),nodci. This is produced by the underestimation of opacities that 
affect the correction for this quantity (Equation [15]). The difference between A^(CO)cmission and the expected CO 
column density is about 20%. 

In the following we use the relation between Af(CO)modci and Af(CO)emission to derive a correction to the observed 
CO column densities. We notice that the difference between these quantities does not show a strong dependence in 
the cloud surface temperature. This is because all models have the same temperature at the cloud center. Since the 
observed excitation temperatures lie between '^ 9— 15K for Ay > 10 mag (Figure [23]). and the excitation temperature 
derived from ^^CO is similar to the kinetic temperature at the cloud surface, we average all models from Tg = 9K to 
15 K in steps of 1 K to derive a correction function. In Figure[24]we show the difference between the model and derived 
CO column density (A^(CO)diff = .^^(CO)modci ~ -^(CO)emission) as a function of the model CO column density. To 
this relation we fit a polynomial function given by 
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Figure 23. {upper row) Excitation temperature derived from the ^^CO emission as function of N{CO)inodal ^or isothermal cloud models 
(left ) and models with temperature gradients (right). In both panels we also show the model cloud kinetic temperature averaged along 
the line of sight, (T]jin)ios, as a function of Af(CO)niodol- (middle row) The -'^''CO opacity versus A''(CO)niodcl for the models shown in the 
upper row. The opacity was derived using Equation II20I I. We also show the opacity calculated from the model cloud N(^^CO) with the 
assumption of LTE. (lower row) Column density of CO calculated from the ^^CO and ^"^CO emission using Equation I I17I I versus the model 
cloud N{CO). The straight line corresponds to an one-to-one relation. 
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Figure 24. Difference between the expected Af(CO)niodcl and the derived A'(CO)emission as a function of the model CO column density. 
The red line represents a second-order polynomial fit. 

To apply this correction we made a roug h estimat e of th e gas-phase CO column density as a function of visual 
extinction. We use the observations bv Frerking et all ()1982D of C^^O and the rarer isotopic species C^^O and ^"^C^^O 
in the direction of field stars located behind Taurus. We convert the observed column densities into A^(CO) assuming 
[CO]/[Ci80]=557, [Ci80]/[Ci^O]=3.6, and [Ci«0]/[i3Ci80]=69 (|Wilsonlll999n . Due to their low abundances, these 
species are likely not affect ed by saturation. No te that, however, they are still sensitive to the determination of 
the excitation temperature. iFrerking et al.l ()1982l ) presented column densities as lower limits when ^^CO is used to 
determine Tox("^^CO) (average ~10K) and as upper limits when they used Tox(^^CO)/2 (i.e. ^ 5K) as the excit ation 
temperature. The kinetic temperature in dense regions is likely to be in between 5 and 10 K (jGoldsmithl 1200 It) and 
therefore, assuming that the isotopologues are thermalized, the excitation temperature should also have a value in 
this range. Thus, we use the average value between upper and lower lim its of the CO column density to determine its 
relation with Ay. For the visual extinction at the positions observed bv IFrerking et al] (|1982f ). we use updated values 
derived bv lShen ov et a l.l (120081) fron i infrared observationE0. We note that the visual extinction correspond to a single 
star while the IFrerking et aLl i 1982D observations are averaged over a 96" beam. W e constructed an ext inction map 
of Taurus with 96" resolution and an extinction curve that matches that adopted by iShenov et al.l (|2008l ) in order to 
compare with their determination of Ay. We found that the visual extinctions always agree within ±0.4 mag. 

In Figure [25l we show the relation between N{CO) and Ay with and without the correction for the effects of temper- 
at ure gradients al o ng th e line-of-sight. For reference we include the values of iV(CO) derived from the observations 
bv IFrerking et all ()1982[ ). The error bars denote the upper and lower limits to the CO column density mentioned 
above. Although our determination of the gas-phase A^(CO)/Av relation is necessarily approximate, the validity of 
the correction for the effects of temperature gradients along the line-of-sight is confirmed by the good agreement 
between Ay and 7V(C0) up to Ay ~23mag after the addition of the column density of CO-ices (Section lS.l.ip . 
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gradients. 
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